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Abstract. We study electromagnetic streaming instabilities in thermal viscous regions of 
rotating astrophysical objects, such as, magnetized accretion disks, molecular clouds, their 
cores, and elephant trunks. The obtained results can also be applied to any regions of in- 
terstellar medium, where different equilibrium velocities between charged species can arise. 
We consider a weakly ionized multicomponent plasma consisting of neutrals and magnetized 
electrons, ions, and dust grains. The effect of perturbation of collisional frequencies due 
to density perturbations of species is taken into account. We obtain general expressions for 
perturbed velocities of species involving the thermal pressure and viscosity in the case in 
which perturbations propagate perpendicular to the background magnetic field. The disper- 
sion relation is derived and investigated for axisymmetric perturbations. New compressible 
instabilities generated due to different equilibrium velocities of different charged species are 
found in the cold and thermal limits when the viscosity of neutrals can be neglected or is 
important. The viscosity of magnetized charged species is negligible for considered perturba- 
tions having wavelengths much larger than the Larmor radius of species. At the same time, 
the neutrals are shown to be immobile in electromagnetic perturbations when their viscosity 
is sufficiently large. 
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1. Introduction 

In a series of papers, Nekrasov (2007, 2008 a,b, and 2009 a,b), a general theory for 
electromagnetic compressible streaming instabilities in multicomponent rotating magnetized 
objects such as accretion disks and molecular clouds has been developed. For accretion 
disks, the different equilibrium velocities of different species (electrons, ions, dust grains, and 
neutrals) have been found from the momentum equations taking into account anisotropic 
thermal pressure and collisions of charged species with neutrals (Nekrasov 2008 a,b, 2009 b). 
For molecular clouds, their cores, and elephant trunks, the different equilibrium velocities 
of charged species have been supposed due to the action of the background magnetic field 
(Nekrasov 2009 a). It has been shown that differences of equilibrium velocities generate 
compressible streaming instabilities having growth rates much larger than the rotation fre- 
quencies. New fast instabilities found in these papers have been suggested to be a source of 
electromagnetic turbulence in accretion disks and molecular clouds. 

Streaming instabilities can also be generated in neutral media. Recently, the streaming 
instability originated due to the difference between equilibrium velocities of small neutral 
solids and gas has been suggested as a possible source that could contribute to the plan- 
etesimal formation (Youdin & Goodman 2005). These authors have numerically treated this 
hydrodynamic instability in the Keplerian disk for the interpenetrating streams coupled via 
drag forces. The growth rates of instability have been found to be much smaller than the 
dynamical timescales. The particle density perturbations generated by this instability could 
seed the planetesimal formation without self-gravity. 

In papers by Nekrasov devoted to electromagnetic streaming instabilities, the viscosity 
has not been considered. However, numerical simulations of the magnetorotational instability 
show that this effect can influence the magnitude of the saturated amplitudes of perturba- 
tions and, correspondingly, the turbulent transport of the angular momentum (e.g., Pessah 
& Chan 2008; Masada & Sano 2008). Interstellar medium, molecular clouds, and accretion 
disks around young stars are weakly ionized objects, where collisional effects play the dom- 
inant role. The ratio of the viscosity to the resistivity (the magnetic Prandtl number) for 
astrophysical objects takes a wide range of values. In particular, in accretion disks around 
compact X-ray sources and active galactic nuclei, the magnetic Prandtl number varies by 
several orders of magnitude across the entire disk (Balbus & Henry 2008). The viscosity of 
astrophysical objects is studied mainly in the framework of the one-fluid magnetohydrody- 
namics. Therefore, one needs to consider this effect using a multicomponent fluid approach 
for studying streaming instabilities in real multicomponent media. 

In the present paper, we investigate electromagnetic streaming instabilities in rotating 
astrophysical objects, such as, magnetized accretion disks, molecular clouds, their cores, 
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elephant trunks, and so on, taking into account effects of collisions, thermal pressure and 
viscosity. We consider a weakly ionized multicomponent plasma consisting of electrons, ions, 
dust grains, and neutrals. The charged species are supposed to be magnetized, i.e., their 
cyclotron frequencies are larger than their orbiting frequencies and collisional frequencies 
with neutrals. We take into account the effect of perturbation of collisional frequencies 
due to density perturbations of species, which takes place at different background veloci- 
ties of species. We find the general expressions for the perturbed velocities of species for 
perturbations which are perpendicular to the background magnetic field. (The analogous 
investigation of perturbations along the magnetic field has been performed in (Nekrasov 2009 
c).) These expressions for any species also contain the perturbed velocity of other species due 
to collisions. The dispersion relation is derived and investigated for axisymmetric perturba- 
tions. The cold and thermal regimes are considered when the pressure force is negligible or 
dominates the inertia. The role of the viscosity of neutrals and charged species is analyzed. 
The growth rates due to different azimuthal velocities of charged species are found. 

The paper is organized as follows. In Section 2 the basic equations are given. In Section 
3 we shortly discuss the equilibrium state. General solutions for the perturbed velocities of 
species for the transverse perturbations having the radial as well as azimuthal wave number 
are obtained in Section 4. These expressions for axisymmetric (radial) perturbations and 
magnetized charged species are given in Section 5. In Section 6 the dispersion relation is 
derived and its solutions are found. Discussion of the obtained results and their applicability 
are given in Section 7. The main points of the paper are summarized in Section 8. 



2. Basic equations 

In the present paper, we investigate weakly ionized astrophysical objects or their def- 
inite regions, where all the charged species are magnetized, i.e., when, in particular, their 
cyclotron frequencies are larger than their collisional frequencies with neutrals. Then the 
momentum equations for charged species and neutrals including the viscous terms have the 
form (Braginskii 1965), 

^ + v r Vv 3i = - W - ^ + ^E ± +^Mv, xB) ± -v jn (v j± -v n± ) (1) 
dt nijiij rrij rrijC - L 

+/XjlVV • Vj + /ij2V 2 VjXZ+ / U 7 - 3 V 2 V ? _ L , 

' h ,z + v, • Vv jz = ^E z +-^- (vj x B) z -v jn (v jz - v nz ) + Afi j3 V\ g , (2) 



at rrij rrijC 

dt 



+ V n • Vv n = -VU - - V U nj (V n - Vj) + fl n ( V 2 v„ + iw • v r 

m n n n . \ 6 



(3) 
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where the index j = e,i,d denotes the electrons, ions, and dust grains, respectively, and the 
index n denotes the neutrals. In Equations (l)-(3), qj and Tfij ^i £1X6 the charge and mass of 
species j and neutrals, Vj >n is the hydrodynamic velocity, n^ n is the number density, Pj >n = 
n j,nTj, n is the thermal pressure, Tj n is the temperature, Vj n = , jj n m n n n (u n j = / ~fj n mjnj) 
is the collisional frequency of species j (neutrals) with neutrals (species j), where 7j„ =< 
av >j n j (mj + m n ) (< av >j n is the rate coefficient for momentum transfer). The coefficients 
of the kinematic viscosity are /iji = Cji/ij/3, = Cj 2 ^j/u} C jTj n , and /ij 3 = Cj 3 /ij / 'u 2 jT 2 n . Here 
f-ij, n — v Tj,nl v 'jn,rm {y nn is the neutral-neutral collisional frequency), u C j = qjB Q /mjC is the 
cyclotron frequency, Tj n = vj^, and vrj, n = (Tj^/rrij^) 1 / 2 is the thermal velocity. Numerical 
coefficients for the electrons and ions are c ei = 0.73, en = 0.96, Cj 2 = 0.5, c e3 = 0.51, and 
Q3 = 0.3 (Braginskii 1965). Further, U = —GM/R is the gravitational potential of the 
central object having mass M (when it presents), R = (r 2 + z 2 ) 1 ^ 2 , G is the gravitational 
constant, E and B are the electric and magnetic fields, and c is the speed of light in vacuum. 
The magnetic field B includes the external magnetic field B$ ext of the central object and/or 
interstellar medium, the magnetic field B$ cur of the stationary current in a steady state, and 
the perturbed magnetic field. We use the cylindrical coordinate system (r,0,z), where r is 
the distance from the symmetry axis z, and 9 is the angle in the azimuthal direction. The 
background magnetic field is assumed to be directed along the z axis, B = B 0zext + B 0zcur , 
the sign _L denotes the direction across the z axis. In Equations (1) and (2), the condition 
oj C j Vj n is satisfied in the viscous terms. For unmagnetized charged particles of species 
j, u C j Uj n , the viscosity coefficient has the same form as that for neutrals. We adopt the 
adiabatic model for the temperature evolution when Pj n ~ n ]% where 7 a is the adiabatic 
constant. In general, we consider the two-dimensional case in which V = (d/dr,d/rd6,0). 

The other basic equations are the continuity equation, 

Oti • 

_p + V . n . nVjn = o, (4) 

Faraday's equation, 



and Ampere's law, 



V x B = —j, (6) 

c 

where j = Y2j Qj n j v j- We consider the wave processes with typical time scales much larger 
than the time the light spends to cover the wavelength of perturbations. In this case, one 
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can neglect the displacement current in Equation (6) that results in quasineutrality both for 
the electromagnetic and purely electrostatic perturbations. 



3. Equilibrium 

We suppose that in equilibrium the electrons, ions, dust grains, and neutrals rotate in 
the azimuthal direction of the astrophysical object (accretion disk, molecular cloud, its cores, 
elephant trunk, and so on) with different, in general, velocities i>j, n o- These velocities can 
depend on the radial coordinate. The stationary dynamics of charged species is undergone 
by the action of the background magnetic field and collisions with neutrals. In their turn, 
the neutrals also experience a collisional coupling with charged species influencing their 
equilibrium velocity. Some specific cases of equilibrium have been investigated in papers by 
Nekrasov (2007, 2008 a,b, 2009 b). Stationary velocities of charged species in the background 
electric field and gravitational field of the central mass in the absence of collisions have been 
used in (Nekrasov 2007). The case of weak collisional coupling of neutrals with light charged 
species (electrons and ions) and weak and strong collisional coupling between neutrals and 
heavy dust grains has been considered in (Nekrasov 2008 a). Equilibrium where the neutrals 
have a strong coupling with light charged species as well as with heavy dust grains has been 
investigated in (Nekrasov 2008 b and 2009 b). 

In papers cited above, it has been shown that the different charged species have different 
stationary velocities. Due to this effect, the electric currents exist in the equilibrium state, 
which generate their magnetic fields. 



4. Linear approximation: General expressions 

In the present paper, we do not treat electromagnetic perturbations connected with the 
background pressure gradients. Thus, we exclude the drift waves from our consideration. 
We take into account the induced reaction of neutrals on the perturbed motion of charged 
species. The neutrals can be involved in electromagnetic perturbations, if the ionization 
degree of medium is sufficiently high. We also include the effect of perturbation on the 
collisional frequencies due to density perturbations of charged species and neutrals. This 
effect emerges when there are different background velocities of species. Then the momentum 
equations (l)-(3) in the linear approximation take the form, 

^+v jr Vv jU = -c^ + ^Eu+^^xB^+^^xBo), (7) 
dt 3 n j0 i ii j in jC x in :i c 
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-v]n (Vjl± - V„ij.) -v]^ (Vjo - V n0 ) + /ijlVV • Vji 



™n0 



+/i j 2V 2 v i iXz+ / u i3 V 2 v j i_ L , 



dvjiz 
dt 



+ v io • Vuj-u = -^£i 2 +^- (vj-o x Bi) -i/? - u nU ) + 4^- 3 V 2 Uju, (8) 



rrijC 



dv. 



nl 



dt 



+ v„ • Vv„i = 



n n o 



n j0 



(9) 



+fin ( W nl + -VV ' V„i 



where c S j jn = (7 a 7j>o/' m j,n) 1//2 is the sound velocity (7 a = 2(3) for the two(one)-dimensional 
perturbations), Uj n = •yj n m n n n0 , zA = ^j n mjnj . The terms proportional to n n \ in Equation 
(7) and riji in Equation (9) describe the effect of perturbation of the collisional frequencies 
u }n = v{ jn{ n ni/n n o) and v^j = ^nj( n ji/ n jo) due to number density perturbations. The 
index 1 denotes quantities of the first order of magnitude. The neutrals participate in the 
electromagnetic dynamics only due to collisional coupling with the charged species. On 
the left hand-sides of Equations (7)- (9), we do not take into account the terms of the form 
Vj,nvVvj tTI Q. These terms without neutral dynamics have been involved in papers by Nekrasov 
(2007 and 2008 a). The neutral dynamics has been included in papers (Nekrasov 2008 b 
and 2009 b). In all specific cases that were considered the terms mentioned above were 
negligible. From general expressions for perturbed velocities given in (Nekrasov (2008 b and 
2009 b) it can be seen that these Coriolis terms can be neglected under sufficient condition 
^j,nz 2Qj,n, where oj^ nz in the present case are given below and Qj jTl are the rotation 
frequencies of species in equilibrium. 

The continuity Equation (4) in the linear regime is the following: 

^wr~ + Vj,n0 ■ Vn 3> i + n iin0 V • v j>n i = 0. (10) 



We further apply the Fourier transform to Equations (7)-(10), supposing perturbations 
of the form exp(ik r r + im6 — icut). Below, we find expressions for the perturbed velocities of 
neutrals and charged species. 
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4.1. Perturbed velocity and number density of neutrals 

Using the z component of Equation (9), we find the Fourier amplitude of the induced 
longitudinal velocity of neutrals i> n izfc, 



iu; nz v, 



nlzk — ^] VnjVjlzk, (H) 



where 

u nz = uj n + iv n + ifi n k 2 j_, 

L0 n = LU - k • V n0 , 

j 

k = {k, u} = {k r , kg = m/r,u}, k\ = k% + kj. Here and below, we omit, for simplicity, the 
index for v® n and 

From Equations (9) and (10), we obtain the induced transverse velocity and induced 
perturbed number density of neutrals, 

- iUnzVnll-k = Vnj^ 2 ^ kr) n — - (v n0 - V j0 ) + ^nj^jlLk, (12) 



n i0 



^nlfc V"^ n-jlk 



J 



tUJ n± = f„j— (13) 

n n0 ^ n j0 
where 

We see from Equation (12) that the perturbed velocity of neutrals is induced by the per- 
turbed number density and velocity of charged species. Equation (13) shows that the density 
perturbation of neutrals is defined by the density perturbations of charged species. When 
u n ± ~ v n j, we obtain n nlfc /n n0 ~ ^jifc/^jo- Note that to obtain the correct expressions (12) 
and (13) it is necessary to take into account perturbation on the collisional frequency of 
neutrals with charged species. 
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4.2. The longitudinal perturbed velocity of charged species 

Let us now find the longitudinal perturbed velocity of charged species Vji z . From Equa- 
tions (5) and (8) in the linear approximation and by using Equation (11) we obtain, 

UJ nz UJj z Vj lzk + Vj n ^2 v nlVl\zk = iu nz F jlzk , (14) 

i 

where 

u jz = LOj + iv jn + i4/i j3 k 2 ± , 
ujj = uj — k • Vjo, 

-Tjlzk — &\zk- 

rrij to 

Equation (14) shows that due to collisions of neutrals with charged species and charged 
species with neutrals, the charged species experience the feedback on their perturbations. 
If neutrals collide with two charged species only, for example, with electrons and ions, the 
solutions of Equation (14) for e and i have the form, 

D z V e \ zk = ia iz F elzk — iVenV n iFii zk , (15) 
D z Vi\ zk iOl ez Fi\ zk ih 'i n h ' ne F elzfc ; 

where 

Oij z U! nz U!j z -\- fj-nMnit 
D z — U nz (jJ iz UJ ez + U iz V en V ne 

^ez^ivMni- 

Substituting Fji zk in Equations (15), we obtain, 

D z Vjizk = i (ajz ~ n g b jz ) E lzk , (16) 

where rig = kgc/u and 



&ez 


Q-iz 

m e 


Qi 

^en ^ni 

rrii 


i 




bez 


Qe 

®-iz 

m e 


VeO 
C 


qi_ 

rrii 


ViO 

c 


Ojiz 


Qi 


m e 


i 




h z 


ft 

O^ez 

rrii 


ViO 

^in^ne 

C 


qe_ 

m e 


VeO 
C 
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If neutrals collide with ions and dust grains, the index e in Equations (15) and (16) 
must be substituted by the index d. 



4.3. The transverse perturbed velocity and perturbed number density of 

charged species 

We now find the transverse perturbed velocity of charged species v^. From Equation 
(7), we obtain two equations for components Vj\ r k and Vjwk, 



lLOj±Vji r k 

-iujj±v jiek 



-lk r C sj \-U cj Vji0 k + (sjlrk 

n j0 



n j0 



(17) 



Here, the following notations are introduced: 



UJ j± 

* 

*2 

Gjlrfik 
Fjlrk 

FjWk 

Qjlrk 
QjWk 



UJ, 



c sj - lUjfiju 

Fjlrfik + Qjlrfiki 



E\ r k-\ — — {n r E iek 



n e E 



Irk) 



rrij 



VjnVnWk — Vjn [VjO ~ V n0 ) . 

n n0 



In the expression for Fji r k, we have used Equation (5). Solutions of Equations (17) have the 
form, 



Dj±Vj lrk = iujj± (1 — Sjog) Gji r k + i (ico^ + ujj±5j r g) Gjwt, 
DjiVjWk = iiOj± (1 — 5j rr ) Gjigk + i (— + UJ j±^jre) Gji r k, 



(18) 



where 



Dj± = ujI (1 - S j± ) - 



CO 



*2 
cj 1 
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and Sjim = kik m c*'j/ujUj±, l,m — r,9, 5j± = 5j rr + 5jgg. The number density perturbation, 
n jik = n jok " Vjik/wj, is the following: 

UjDj_J^ = (ik r u j± + k e u* C] ) G jlrk + (ik g uj j± - k r u* CJ ) G jWk . (19) 

In the absence of the viscosity, expressions (18) and (19) coincide with the corresponding 
expressions in the paper by Nekrasov (2009 b). 



5. Axisymmetric perturbations 



Below, we shall consider axisymmetric perturbations with k r ^ 0, kg = 0. In this case, 
Equations (18) and (19) take the form, 



D j± v 



j± Ujlrk 



Qj 



= iuOj±—E lrk H — — ( iuj 



Qj 



k r VjO 



CO 



U*j ) E Wk + iuJjxUjnVnirk 



(20) 



n n \k / \ 

VnWk [VjO ~ V n0 ) 



D j±Vj 



19k 



ft x \ , * krV i° 

lL)j± (1 — Ojrr) + L) cj 



?J p . * Qj 771 . * 

&\g k + 0J c j hii rk + OJ c jVj n V n lrk 



m 



rrii 



-\-%UJj\_Vj n (1 — $jrr) 



^nlk I \ 

VnWk (VjO ~ V n0 ) 

TlnO 



Tlj\ k k r Vj\ rk 



n j0 



(21) 



From Equation (12), we can find components of the perturbed velocity of neutrals in the 
present case, 



VnXrk 



1 / u VnjVjlrk-, 



(22) 



VnWk — 



i — y v nj — [v n0 - Vjo) + % — y 



VnjVjWk- 



Substituting these expressions in Equations (20) and using Equations (13) and (21), we 
derive equations that contain only perturbed velocities of charged species, 



Dj±v 



j-L Ujlrk 



Qj 



= iuiji — —E\ rk H — — iuj 



Qj 



J j±- 



k r v j0 



u*j I E wk 



Mirk 



(23) 
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D j± Vj 



16k 



—I 



00 nz 



+1 



tUJj± (1 — 5j rr ) + W CJ 



V in U) 
CO 



&r (V n0 -Vio) k r (v j0 ~ V n0 ) 



+ 



00 r . 



00 n ± 



00 n ±_ 



E 



LO 

OJj±Vjn 



1l 771 ,,,* ?3 p 



(1 - £ J>r ) ^ Vn, 



iVlWk 



OOj\_Vj n 
00 



1 $jrr) ^ ^ ^Vi 



iVllrk 



K (V n0 - V i0 ) . K (VjO ~ V n0 ) 



Mr. 



+ 



0->n± 



We shall find solutions of Equations (23) for the case in which charged species are 
magnetized, i.e., when the Lorentz force dominates inertia, thermal, and drag forces. More 
specifically, we suppose that 



00j± 



LO 



(1 — 5j rr ) <C 



KVjO 



00. 



LO 



LOn 



(24) 



'1 — Sj rr ) must be satisfied. 



It follows from these inequalities that the condition lo^ oo? A 
From the latter, we obtain that oo^j oo^ ± and oo/oOj± k^Pj, where pj = c S j/LO C j is the 
Larmor radius. Thus, the viscosity of charged species will be negligible for our consideration 
given below. Besides, c*| ~ c 2 sj because we consider that lo <C u jn . In the case (24), the main 
perturbed velocities of charged species (denoted by the upper index 0) are the electric and 
magnetic (due to the Lorentz force) drifts, 



u jlrk 



J jWk 



rrijLOcj 



-E 



iek, 



(25) 



rrijLO C j 



k r Vjo 
£>irk H £>Wk 



LO 



Substituting expressions (25) in the collisional terms in Equations (23), we obtain the fol- 
lowing expressions for the transverse velocities of charged species: 



Dj±Vji r k 

Dj±VjWk 



m, 



LO 



LO 



. <Ij <lj I KVjO . . k r V n o\ . ■ 
l — CTjrr&lrk + — [ —<7 jrdl + 1>Gjr62 ~ r ?0>6>3 ] &16k 

rrii 



GjQ r E lrk 



rrti 



1>&jr62 ~ ttOj±Oj rr + \<Jj r ei ~ & jOBl) r (TjQQ 2 

LO LO 



Ewk, 



(26) 



where notations are introduced, 



' jrr 



= 00 i± + 



OOnz 



00jl_VjnVn 

, &jr61 — 00 C j H , <Jjr62 = 00 j± + 

00 n ±00 c j 



OJ n ± 
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O"jr>03 — ^jn^n I I j °j0r — ^cj H U _ Ojrr) , 

UJj\_Vj n V n _ Uji_l>j n l> n . ( 1 1 

Cj001 — - 



-Oj rr , <Xj002 — \ Y — Ojrr) I I 

LO c j \ UJ nz UJ n j_ J 



Below, we shall derive the dispersion relation and find its solutions. 



6. Dispersion relation 

From Faraday's and Ampere's laws (5) and (6), we find the following equations: 



J Irk 



0, 

Ani 



(27) 



n^Eiok = ji0 k , 



n 2 r E lzk 



uj 
Ani 

uj 



-Jlzki 



where j lr , zk = Y,j qjnjov jlr , gk , jwk = J2j QjrijoVjiek + J2j QjnjikVjo, and n 2 = k 2 r c 2 /uj 2 . We see 
that perturbations with polarizations of the electric field along and across the background 
magnetic field are split. We do not consider perturbations with longitudinal polarization 
because the current ji zk does not contain the equilibrium velocities in the axisymmetric case 
(see Equation (16)). Using expressions (26), we find the transverse electric currents, 



Ani . 

UJ 

Am . 

Jwk = £e r Eirk + ZeeEwk- 

UJ 



(28) 



Here, 



£eo — 



ST rr r 

3 J 



E 

3 

E 



UJ 



pj 



u)Dj\ 
u)D,j\ 



E 



UJ 



PJ 



ujD 



3± 



KJjQr ~ °" 



k r Vjo 



jrr 



UJ 



k r Vjo k r v n o 

KTjrOl + <Jjr62 r 0>6>3 



UJ 



UJ 



i i i kr v jo\ r ( ■ k r Vj \ k r v n0 

-<Jjr62 \ 1 H — — I + UJj±0jrr + I M? jee2 ~ ^jr93' 



UJ' 



UJ 



UJ 



where uj p j 



{AitrijQq 2 /mj} 1 ^ 2 is the plasma frequency. 



-13- 



From Equations (27) and (28), we obtain the dispersion relation, 

e rr n 2 = E rr Sgg — E r gEg r . (29) 

Below, we will find solutions of Equation (29) in limiting cases. 

6.1. The case of cold species 

We at first consider the case in which neutrals and charged species are sufficiently cold, 
uj 2 ^> k 2 c 2 n and uvj n ^> k 2 c 2 j (it is obvious that the first condition is the main one). We 
also suppose that Uj n ^> v n ^> uj. Under these conditions the values a have a simple form, 

Oj„ = CTj r 02 = ^Vjn/lSn, O jg r = <J j r 0i = UJ c j + V^J 'uJ c j + iuJV 2 r J LU c jV n , and (Jj r 03 = <Jjgg 2 = 0. 

When calculating the value e rr Egg — e r geg ri we have carried out symmetrization, used the 
condition of quasineutrality, ^ ■ qjUjo = 0, and kept the main terms. Then the dispersion 
relation (29) takes the form, 



U ^ -J^-J- ^J~ k r C ~ 2 W 2. ^2 *r (^0 " «») ■ (30) 

j7 C J cl j ^cj * j t ^c] ^cl 

We see that the streaming instability is possible when the difference of the equilibrium 
velocities of charged species is sufficiently large to exceed the threshold defined by the first 
term on the right-hand side of Equation (30). 

Let us consider, for example, Equation (30) for the electron-ion plasma (without dust 
grains). Then this equation will be the following: 

2 ; 2 2 m e l 'en 1 2/ \2 /oi\ 

LO = k r C A k r {V e0 ~ Via) , (31) 

1 /9 

where ca = (B 2 / A7rm n n n0 ) is the Alfven velocity. When obtaining this equation, we have 
used that miU in ^> m e u en (see e.g. Section 7). The first term on the right hand-side of 
Equation (31) describes the magnetosonic waves in a weakly ionized plasma. The second 
term can generate the streaming instability of the hydrodynamic kind. This term does not 
change the order of the dispersion relation for perturbations under consideration that results 
in appearance of the threshold for the absolute instability to develop (see also Nekrasov 
2007). 
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6.2. The case of thermal species 

In the case of thermal species, we suppose conditions fi n k 2 ^> uj and kfc 2 j ^> ujVj n to be 
satisfied. Then we have that k 2 c 2 n ^> uj 2 (the intermediate case for neutrals, fi n k 2 <C uj <C 
k r c sn , we will not consider). We also adopt, as above, that vj n u n . However, the condition 
u n ^> uj is now not necessary. Then the values a take the form, 



IKV 



fink 2 . 



jn 



(Jjrdl — W c j, CTj r 02 — IV jn-, 0>6»3 — — 1>KVj n , 



where k = v n j {v n + fi n k 2 ) and dj = k 2 c 2 j/ujuj C j. Substituting these values in e rr ^ and eg r ^, 
we find expression e rT eee — £ r e£er, 

fi n k 2 



£rr£06 ~ £r6£0r — lK / AjAjGtj + K £ AjAj 



jl 



k 2 . (Vj - vmY 
2lu 2 



(32) 



where Xj = —uj^-Vj n jujuj 2 c -. When deriving Equation (32), we have used that ^ . oj 2 jUJ C j/ujDjj_ = 



—i ^ . Xjdj due to condition of quasineutrality. The first term on the right hand-side of this 
equation is the main one according to conditions at hand among other terms that are inde- 
pendent of background velocities. We see from Equation (32) that the viscosity of neutrals 
plays an important role in the case under consideration. 

The dispersion relation (29) takes the form, 



UJ 



2 £ 



2 



UJ 



ji 



2 2 



{vinC% + V jn C 2 sl ) 



—2 77T- ( v jo - v io) 



jl 



UJ 



UJ 



(33) 



cl 



Thus, the streaming instability exists for the sufficiently short wavelengths of perturbations 
when the viscosity of neutrals and thermal pressure of species play a role. 

We now consider Equation (33) for the electron- ion plasma (as above). In this case, we 
obtain, 

.m e (v e0 - v i0 ) 2 



UJ 



m i i C Ai + C s) 



(34) 



where ca% = {Bl/ -iTrniiTiio) 1 ^ is the ion Alfven velocity and c s = [y a (T e + T,j) /rrii] 1 ^ 2 is 
the ion sound velocity. We see from Equation (34) that the growth rate 7 =Im uj can be 
sufficiently large in comparison to the rotation frequency in media where c 2 Ai + c 2 is not too 
large as compared to (v e o — Vio) 2 . 
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The instability described by Equation (33) belongs to the class of dissipative instabilities 
because the dispersion relation contains the corresponding imaginary term which is propor- 
tional to the collisional frequency Without additional damping mechanisms, the threshold 
of this dissipative-streaming instability is absent. 

7. Discussion 

Let us discuss some consequences that follow from the results obtained above. From 
Equation (13), we see that the relative number density perturbation of neutrals is of the 
order of that for charged species, if u nl _ ~ u n . However, in the case k 2 c 2 sn ^> u)v ni a neutral 
fluid can be considered as incompressible. In the latter case, the perturbed radial velocity 
of neutrals in the axisymmetric perturbations is also negligible (see Equation (22)). At the 
same time, the perturbed azimuthal velocity of neutrals is of the order of that for charged 
species for small viscosity of neutrals, v n ^> p n k 2 {y n ^> u), and is negligible for large 
viscosity of neutrals, v n <C p n k 2 (see Equation (22)). Thus, the neutrals are immobile in 
perturbations when their viscosity is large. 

As a specific example, we consider the case of the dense regions of molecular cloud 
cores where the number density of neutrals is n n ~ 10 5 cm~ 3 and rii/n n ~ 1CT 7 (Caselli 
et al. 1998; Ruffle et al. 1998). For simplicity, we omit the index by rij . The ratio 
of mass densities of dust grains, pd = rridrid, and neutrals, p n = m n n n , in the interstellar 
medium is of the order of 1CT 2 (e.g. Abergel et al. 2005). From these relations, we obtain 
that rid <S rii at m d 10 5 m n . For usually adopted m n = 2.33m p , where m p is the proton 
mass, and for ^ = 1 g cm~ 3 , where is the density of grain material, we obtain that 
rrid ^> 10 5 m„ = 3.9 x 10~ 19 g is satisfied at ^> 4.53 x 10~ 3 /im, where is the grain 
radius. The typical values of grain radius are ~ 0.01 — 1 pm (Mendis & Rosenberg 1994; 
Wardle & Ng 1999). However, in spite of that rid <C rii the mass density of dust grains 
p d ~ 7.77 x 10 3 pi for m, = 30m p . 

Further, we take the magnetic field B = 150 pG and the radius and charge of dust 
grains = 0.01 pm and qd = q e (we consider one typical type of dust grains). In this 
case, rrid = 4.19 x 10~ 18 g and n d = 0.93 x 10~ 3 cm~ 3 . Then we obtain the following 
values for the plasma and cyclotron frequencies of charged species: u pe = 5.64 x 10 3 s _1 , 
tj ce (> 0) = 2.64 x 10 3 s -1 , LU P i = 24 s _1 , u ci = 4.78 x 10~ 2 s" 1 ($ = -q e ), u pd = 2.53 x 10~ 2 
s -1 , and u) c d{> 0) = 5.75 x 10 -7 s -1 . The rate coefficients for momentum transfer by elastic 
scattering of electrons and ions with neutrals are < av > en = 4.5 x 10~ 9 (T e /30 K) 1 / 2 cm 3 
s- 1 and < av > in = 1.9 x 10~ 9 cm 3 s" 1 (Draine et al. 1983). We take T e = 300 K. Then we 
obtain u en = 1.42 x 10~ 3 s _1 and u in = 1.37 x 10~ 5 s -1 . The collisional frequency of dust 
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grains with neutrals has the form v dn ~ 6.7 ' p n r 2 d VTn/'md (T n ~ T e ) (e.g., Wardle & Ng 1999). 
Using parameters given above, we obtain u dn = 6.24 x 1CT 8 s _1 . The collisional frequency of 
neutrals with charged species is v n ~ v nd = 6.24 x 1CT 10 s _1 . 

The values u 2 jUj n /uj 2 j which are contained in Equations (30) and (33) are under condi- 
tions at hand the following: uj* e v en /ujl e = 6.4 x 10~ 3 s" 1 , uj^Vin/uj 2 ^ = 3.44 s" 1 , uj 2 d v dn juj 2 cd = 
120.8 s _1 . Thus, we can write Equation (30) in the form, 



2 ^nd 
UJ = 

Vdn 



"in (V i0 ~ V d0 ) 



Vnd 



-Ai 



k 2 r 2 



where CAd = cuo cd /uo pd is the dust Alfven velocity. Supposing that \vio — v d o\ ^> [y n dl '^inY^CAi 
(the sign | | denotes an absolute value), the growth rate of instability 7 will be equal to 



7 



\ 1/2 

Vin \ CAd , I 

k r \v i0 

Vdn J CAi 



v d o\ ■ 



(35) 



This solution exists if 



Vn I V dn \ 1/2 C A i 1 



CAd 



\ 1/2 
, V dn \ CAi 

^dol > C sn \ . 

Vin J CAd 



From these inequalities, it follows that A r ^> 2nc sn /v n) where A r = In jk r is the wavelength of 
perturbations. For parameters given above, we obtain c sn = 1.73 km s _1 ( / ~f a = 3), CAi = 600 
km s -1 , CAd = 6.8 km s _1 , \v i0 — v d0 \ > 10 km s _1 , and A r 1.73 x 10 10 km. The estimation 
of the growth rate (35) at \v i0 — v d0 \ = 15 km s _1 and A r = 5 x 10 10 km is the following: 
7 = 3 x 10~ 10 s _1 . The condition of magnetization (24) is satisfied for ions as well as for 
dust grains. The case of unmagnetized dust grains has been considered in (Nekrasov 2009 
a) 

In the case of thermal species and/or short wavelength perturbations, \i n k 2 ^> u and 
k 2 c 2 j ujUj n , we consider the electron-ion plasma because conditions k 2 c 2 sd ^> u)v dn and 
k r v d o/u <C u)cd/v dn (see inequalities (24)) are incompatible for solution (33) and parameters 
given above. In the absence of dust grains, Equation (33) has the form (34). The condition 
k 2 c 2 si ^> 7z/ in and the right inequality (24) for ions Kvio/j <C oj C i/v in (for electrons these 
conditions are weaker) are compatible if 

mvin (3 uj 2 ci 
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where (5 = Yl-nrii (T e + Tj) /Bq. For number densities rii and n n used above, this inequality 
is also not satisfied. At (3 <C 1 it can be written in the form, 



For n n ~ 10 4 cm 3 and ni/n n ~ 10 6 , this condition is satisfied. 

We have excluded drift waves in our study. This can be done if the frequency spectra of 
these waves and perturbations under consideration are different. As is known, the frequency 
of drift waves, for example, in the electron-ion plasma, is equal to uj — (T e /miUj ci ni) [k x Vnj] 2 . 
For axisymmetric perturbations, this frequency is equal to zero, if the number density is uni- 
form in the azimuthal direction. If dni/riirdO ~ L e 1 ^ in any region, then uj ~ uj C ik r pfL e l . 
Comparing this expression, for example, with the growth rate (34), we see that under con- 
dition 



(f3 <C 1), the growth rate is larger than the drift frequency. The last condition can easily be 
satisfied because pi is much smaller than L e . 



In the present paper, we have studied electromagnetic streaming instabilities in the 
thermal viscous regions of rotating astrophysical objects, such as magnetized accretion disks, 
molecular clouds, their cores, and elephant trunks. However, the results obtained can be 
applied to any regions of interstellar medium, where different background velocities between 
charged species can arise. 

We have considered a weakly ionized multicomponent plasma consisting of electrons, 
ions, dust grains, and neutrals. The cyclotron frequencies of charged species have been 
supposed to be larger than their collisional frequencies with neutrals. The axisymmetric 
perturbations across the background magnetic field have been investigated. We have taken 
into account the effect of perturbation of collisional frequencies due to density perturbations 
of species. New compressible instabilities generated by the different equilibrium velocities of 
species have been found in the cold and thermal limits either when the viscosity of neutrals 
can be neglected or when it is important. For the perturbations considered, the viscosity of 
magnetized charged species is negligible. 

In dense accretion disks, the ions are unmagnetized while the electrons remain magne- 
tized (e.g., Wardle & Ng 1999). In this case, the viscosity of ions can play the same role as 





8. Conclusion 
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that for neutrals. However, our present model does not describe such objects. 

The electromagnetic streaming instabilities studied in the present paper can be a source 
of turbulence in weakly ionized magnetized astrophysical objects in regions where thermal 
and viscous effects can play a role. 

The insightful and constructive comments and suggestions of the anonymous referee are 
gratefully acknowledged. 
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